Zero-field magnetization reversal of two-body Stoner particles with dipolar interaction 
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Nanomagnetism has recently attracted explo- 
sive attention, in particular, because of the enor- 
mous potential applications in information indus- 
try, e.g. new harddisk technology, race-track 
memory [1], and logic devices [2]. Recent tech- 
nological advances [3] allow for the fabrication 
of single-domain magnetic nanoparticles (Stoner 
particles), whose magnetization dynamics have 
been extensively studied, both experimentally 
and theoretically, involving magnetic fields |449|] 
and/or by spin-polarized currents [10- 20]. From 
an industrial point of view, important issues in- 
clude lowering the critical switching field Hc^ and 
achieving short reversal times. Here we predict a 
new technological perspective: He can be dramat- 
ically lowered (including He = 0) by appropriately 
engineering the dipole-dipole interaction (DDI) 
in a system of two synchronized Stoner particles. 
Here, in a modified Stoner- Wohlfarth (SW) limit, 
both of the above goals can be achieved. The ex- 
perimental feasibility of realizing our proposal is 
illustrated on the example of cobalt nanoparti- 
cles. 

The magnetization dynamics of two Stoner parti- 
cles being subject to DDI and an external magnetic 
field is governed by the Landau-Lifshitz-Gilbert (LLG) 
equation0,|2l|,[22|. 
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Here rhi = Mi/Mg is the normalized magnetization vec- 
tor of the i-th particle, {i = 1,2). Ms = \Mi\ is the 
saturation magnetization of either particle, and a is the 
Gilbert damping coefficient. For simplicity, we are as- 
suming the two particles to be completely identical in 
shape, volume, a, and Mg. The unit of time is set to be 
(|7|M5)~^, where 7 is the gyromagnetic ratio, and the 
total effective field hi is given by hi = —dE/drhi^ where 
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i=l,2 



mf ^ -\- fhi ■ h) -\- r][rhi ■rh2 — 3(mi •n)(m2 -n)] 



(2) 

is the total energy per particle volume V in units of 
fioM^ (/io* vacuum permeability). Here both particles 
have their easy axis (EA) along the z-direction, and 



the parameter k summarizes both shape and exchange 
contributions to the magnetic anisotropy. In addition, 
h = H /Ms where i? is a homogeneous and static exter- 
nal field. The parameter r] = -^^ is a geometric factor 
characterizing the DDI with d being the fixed distance 
between the two particles whose direction is described by 
the unit vector n. Here we omit the exchange interaction 
energy between the two particles since it becomes impor- 
tant only at very small particle distances. Moreover, in 
the synchronized magnetic dynamics to be investigated 
below, it only contributes a constant to the energy and 
will therefore not change the physical behavior. 
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FIG. 1: The (normalized) critical switching field he versus the 
DDI strength 77 for PERP and PARA configuration illustrated 
in the insets. Analytical results (solid and dashed lines) are 
compared with numerical findings (circles and squares). The 
system parameters are /c = 0.5 and a — 0.1. 

Let us focus on two typical configurations where 
the connecting unit vector n is either perpendicu- 
lar or parallel to the anisotropy axes, referred to as 
PERP and PARA configuration, respectively (see in- 
sets in Fig. d]). That is, introducing spherical co- 
ordinates rfii = (sin^i cos^i,sin6>i sin^i,cos6>i), n = 
(sin 6>n cos 0n, sin 6>n sin ^n, cos 6>n), we have On = 7r/2 for 
PERP configuration and On = for PARA (or On = tt 
since DDI is invariant under h \-^ —h). Moreover, with- 
out loss of generality, we can let (/)n = by choosing 
the X-axis along the connecting line of the two particles. 
Furthermore, we concentrate on the synchronized mag- 
netic dynamics of the two Stoner particles, where both 
magnetization vectors remain in parallel throughout the 
motion, Oi = O2 = 0^ 4>i = 4>2 = 0- Thus, the two parti- 
cles behave like a single entity, and this two-body Stoner 
particle can be regarded as a computer information bit. 



We have verified by numerical simulations (see discus- 
sion below) that this dynamical regime is stable against 
perturbations. For this synchronized motion mode, the 
nonlinear coupled LLG equations ([1]) read in spherical 
coordinates 

9 ^ a sin Ocj) = —3/^ cos '^ sin On sin (j) , 

aO - smO(j) = hsinO - ksm20 + ^^^^L_^ ^ (3) 

2 uO 

where tp is the angle between m and n, i.e. cos?/^ = 
cos cos On + sin sin On cos (j). Here we have put h = —hz 
(antiparallel) along the easy axis, which is the usual field 
configuration for reversing a magnetic bit. The above 
equations are the starting point of our numerical calcu- 
lations to be discussed below. 

However, in order to analytically explore the SW limit 
for magnetic reversal [7, 23], we assume the external field 
to lie in the plane spanned by the easy axis and the inter- 
particle direction, h = hzZ -\- hxX. The energy takes the 
formE = — 2/ccos^ 0—2hzCosO—2hxSmO-^r]{l—3cos'^ ^fj) 
with ip = — On- The SW limit occurs at the inflection of 
the energy as a function of 6>, i.e. dE/dO = d'^E/dO'^ = 0. 
An elementary calculation translates this condition into 
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(4) 



where the minus (plus) sign corresponds to the 
PERP(PARA) configuration, respectively. Note that in 
the absence of DDI, 77 = 0, the above equation just recov- 
ers the usual SW limit for a single S toner particle. As a 
result, the critical switching field he (applied antiparallel 
to the easy axis) is given by 



h 
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= \2k-3r]l 



PARA 



= 2k^3r]. (5) 



This analytical solution (J5j) is shown by the solid and 
dashed lines in Fig. [TJ 

Moreover, the above findings imply the remarkable ob- 
servation that in the PERP configuration there exists a 
critical DDI strength rjc = 2k/3 such that the critical 
switching field vanishes! {he = 0). 

The zero-field condition is achieved for interparticle dis- 
tances given by 



dc 
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where K = kji^M'^ is the standard anisotropy coeffi- 
cient. Remarkably, the above condition is independent of 
the damping, and its physical contents can be illustrated 
in terms of energy landscape. In the PERP configura- 
tion, the energy in the absence of an external field is 
E = —2r]-\-3{r] — r]c) cos^ 0. Thus, for r] = rjc^ any angle 
represents an equilibrium position such that an arbitrary 
small field is sufficient to move the magnetization along 



the field axis, implying the possibility of zero-field rever- 
sal. Moreover, for r] > r]c^ the zero-field ground state 
is given by ^ = 7r/2, i.e. the synchronized magnetiza- 
tion points along the interparticle axis, while for r] < r]c, 
the ground state magnetization is along the easy axis, 
l9 = 0,7r. 

In order to complement and quantitatively support our 
previous discussion, we have performed numerical sim- 
ulations of the LLG equations (|3]) using the 4th-order 
Runge-Kutta scheme. We consider a range of the DDI 
parameter of < 77 < 1. The case, 77 = corre- 
sponds to the limit d ^ 00, i.e. the two nanopar- 
ticles being infinitely apart. Large rj can be realized 
by fabricating magnetic particles of ellipsoidal shape al- 
lowing for a closer proximity. Throughout the numeri- 
cal results shown here, we use a damping parameter of 
a = 0.1. In Fig. [T]we compare simulation results for the 
critical switching field with the analytical formulae (J5]), 
both findings being in excellent agreement. The criti- 
cal switching field in the PARA configuration is always 
higher than the value without DDI. Thus, only the PERP 
configuration will be useful for possible technology appli- 
cations. 

Let us now turn to the reversal time T^, i.e. the dura- 
tion of the reversal process of the magnetization direction 
changing the angle from to tt. To avoid the metastable 
points 6> = 0,7r, we introduce two small deviations (5^,(5/, 
by defining Oi = 5i (i.e. rriz ~ 1) and Of = tt — Sf (i.e. 
rriz ~ —1), which leads to a finite reversal time in our 
simulations. In the PARA configuration, one can derive 
a closed analytical expression for this quantity. 
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where ho = 2k ^ 3r] and A = (1 + a'^)/a. Note that the 
case T] = also includes the PERP configuration since 
both configurations are indistinguishable here. 

In the upper panel (a) of Fig. [2] we show simula- 
tion results for Tr in the PARA configuration for again 
k = 0.5, Si = Sf = 0.001, and various rj. The simula- 
tion data is well described by the approximate expres- 
sion Tr — — h2~Lh2 valid for small Si = Sf ^ 1. As 
a result, in the PARA configuration the reversal time 
increases with increasing strength of the dipolar interac- 
tion. The sensitivity of the data to the parameters Si 
(where Sf is fixed or vice versa) is illustrated in the bot- 
tom panel (b) of Fig. [2l In accordance with Eq. ([7]), the 
reversal time depends only logarithmically on the quan- 
tity and does not change its order of magnitude while 
Si is changing over several powers of ten. Thus our re- 
sults are not qualitatively affected by our choice of the 
condition Si = Sf = 0.001. 

A similar result regarding the sensitivity to initial con- 
ditions is obtained in the PERP configuration. Here an 
analytical result comparable to Eq. (O does not seem to 
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FIG. 2: PARA configuration: (a) Reversal time Tr versus field 
magnitude h at different DDI strength 77. The simulation data 
(symbols) is fitted according to the approximate expression 



2Xh{-lnSi) 



for 5i 



0.001, and k = 0.5, a = 0.1 



(cf. Fig[T]). (b) Dependence of Tr on the initial condition 5i'. 
Simulation data along with exact analytical results according 
to Eq. dg. The final state is fixed to be 5f = 0.001. 



be achievable. However, the simulation data is shown 
in Fig. [3fb) that Tr again depends only logarithmically 
on 5i and is therefore similarly insensitive to the initial 
condition as in the PARA case. 

The dependence of the reversal time on the dipolar 
parameter 77 shown in Fig. [3{a), however, is strikingly 
different from the PARA configuration: Here Tr clearly 
decreases with increasing 77, especially in the regime of 
low fields, which also strongly favors the two-body bit 
implementation proposed in this letter. 

Thus, around the critical DDI strength rfc we find not 
only nearly a zero critical switching field but also a sub- 
stantially shorter reversal time. This result promises 
attractive future applications for computer technology, 
such as fast read/ write hard disk or magnetic RAM. 
Let us discuss two schematic setups for experimentally 
realizing the zero-field switching mechanism. The first 
scheme A is illustrated in the upper-left panel of Fig. HI 
along with a numerical simulation in panel (a). Here 
we chose again k = 0.5, and the DDI parameter is 



FIG. 3: PERP configuration: (a) Reversal time Tr obtained 
from numerical simulations versus field magnitude h at differ- 
ent DDI strength 77. (b) Tr as a function the initial condition 
5i. All other parameters are the same as Fig. [21 



T] = 0.32 < Tjc = 1/3. Thus, the magnetization in the 
zero-field ground state points along the easy axis (z-axis). 
Implementation steps of scheme A are demonstrated as 
follows: Step 1 - First we choose a casual initial con- 
dition, for instance miz{0) = m2z{0) = 0.5, to mimic 
the relaxation of the system to its ground state under 
zero-field during < t < 500; Step 2 - Then we apply 
a tiny antiparallel field h = 0.03 during 500 < t < 1000 
which drives the magnetization reversal process of the 
two-body Stoner particle along the field direction. The 
reversal time is found to be Tr ~ 198(|7|Ms)~^, and the 
inset shows the reversal process on a smaller scale; Step 
3 - Finally we quench the field and obtain the new stable 
magnetic state. 

The second scheme B is sketched in the upper-right 
panel of Fig. |4] along with a numerical simulation in mid- 
dle panel (b). Here r] = 0.34 > r]c = 1/3 for k = 0.5, 
i.e. the zero-field ground state is along the hard axis 
(x-direction), and a field along the easy axis is per- 
manently required to preserve the magnetization state 
(information) [24] . However, our results show that such a 
field can be very small and actually close to zero. Thus, it 
is not implausible to generate such a field as the Oersted 
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FIG. 4: Upper panels: Two schemes of magnetization rever- 
sals of the two-body Stoner particle (see text). 
Middle panels: mz{t) versus time at (a) 77 = 0.33 < rfc — 1/3, 
(b) Tf — 0.34 > rfc — 1/3. In (a), an antiparallel field h — 0.03 
is applied only during the time interval 500 < t < 1000 with 
a reversal time of Tr ~ 198(|7|Ms)~^; In (b), a parallel field 
h = —0.03 is applied during 500 < t < 1000 and then reversed 
to /i = 0.03, leading to a reversal time as Tr ~ 155(|7|Ms)~^. 
The insets show the reversal process on a smaller scale. 
Bottom panels (c) & (e): Average stable magnetization 
{miz + 77122 )/2 versus the deviation angle 6 = O2 — Oi (in 
degrees) between the initial magnetization directions in the 
two schemes A and B, respectively. Panels (d) & (f): Anal- 
ogous data as a function of the deviation field 6h = h2 — hi 
antiparallel to the easy axis (see text). 



field of tiny switchable currents. The reversal process can 
be implemented as follows: Step 1 - We again choose the 
initial condition miz{^) = Tn2z{0) = 0.5, and the system 
relaxes to its ground state (x-axis) during < t < 500. 
Then we apply a tiny field h = —0.03 (note our defini- 
tion h = —hz) to preserve the magnetic state rriz = 1 
during 500 < t < 1000; Step 2 - After t > 1000 the 
field is reversed to h = 0.03, and the two-body parti- 
cle reverses its magnetization during a reversal time as 
Tr ^ 155(|7|Ms)"^ (see also inset); Step 3 - The reverse 
field is applied permanently to preserve the bit informa- 
tion. 

We now address the stability of the synchronized dy- 
namics of the two-body particle as studied so far. Let us 
first investigate small deviations of the initial magneti- 
zation directions of the two subparticles being otherwise 
still identical. In detail, we fixed the initial direction 
(at t = 500 in reversal scheme A, and at t = 1000 for 
scheme B) of one subparticle to be close to the z-axis, 
i.e. Oi = 0.001, and changed the other initial direction 
to 6>2 = 0.001 + 5. The result is shown in Fig. gjc) and 
(e): For a finite range of deviations, 5 < 7.8° {5 < 4.5°) 
in scheme A(B), the synchronized mode remains stable, 
while for substantially larger J, the average magnetiza- 
tion mostly reaches zero. 

Let us now consider the case where the two particles 
differ slightly in anisotropy Ki, volume Vi^ saturation 
magnetization Mg^i {i = 1,2). Here the zero-field mecha- 
nism still occurs at r]c = {h + fc2)/3 = 3/xoMAM^2TA+k) ' 
However, the effective field experienced by each particle 
will also be slightly different. To numerically check the 
stability under such different external fields, we fixed the 
field on one subparticle to hi = 0.03 and changed the 
field on the other subparticle as /12 = 0.03 + Sh. The nu- 
merical results are given in Fig IH^d) and (f ) for reversal 
schemes A and B, respectively, and demonstrate again a 
finite range of stability against deviations from the case 
of strictly identical particles. 

To give a concrete and practical example of our find- 
ings, let us discuss the case of cobalt (Co) particles. The 
standard data is Ms = 1400kA/m, uniaxial strength 
K = 10^J/m^ a = 0.1[4]. Thus k = K/{fioM^) = 0.04 
such that r]c = 0.027. For two spherical particles with 
radius r, i.e. r] = r^/{3d^), the critical DDI strength is 
reached at dc = 2.3r. The critical switching field with- 
out DDI (SW limit) is Hsw = 2K/{^qMs) = 1400Oe. 
In the presence of DDI, and considering deviations Ad 
from dc, one can express the critical switching field as 
He/ Hsw = 3\Ad\/dc. Thus, in order to drastically 
reduce the switching field by taking advantage of our 
proposal, one has to engineer the interparticle distance 
on a scale of dc = 2.3r which is typically a few hun- 
dred nanometers. In the case of Co, the time unit is 
(|7|M5)~^ = 3.23ps rendering the reversal times in the 
schemes A and B to be Tr ~ 0.64ns and Tr ~ 0.5ns, 



respectively, which are much shorter than in the conven- 
tional setup. 

Another important issue concerns thermal fluctua- 
tions. For T] < Tjc the ground state (^ = 0, tt) is stable if 
the energy barrier in presence of DDI is large compared 
to the thermal energy 'i{/\r])iiQM'^V > ksT {ks: Boltz- 
mann's constant), which translates for T = 400K and 
the above parameters for Co into {Ar])r^ :^ O.lSnm^. 
On the other hand, the energy scale of the applied 
field should also be large compared to thermal effects, 
2jiioMsHV ^ ksT^ which, under the same conditions, 
reduces to H :^ 4700/r^ Oenm^. Thus, considering a 
typical particle radius of r = lOOnm, i.e. dc = 230nm, 
and Hc/Hsw = 0.03, i.e. He = 420e, both conditions 
are easily satisfied, (A77)r^ = 810nm^ ^ O.lSnm^ and 
He = 420e > 0.0047Oe. 

Finally, we would like to remark that our general re- 
sults regarding the infiuence of dipolar interaction on 
the critical switching field of two-body Stoner particles 
are consistent with recent experimental and micromag- 
netic simulation results showing that the coercive field 
for an array of nanowires can be lower than for a single 
nanowire[25, 26]. 
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